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1 Introduction 



Several integrable quantum spin chain models within the range of the algebraic Bethe ansatz 
method have a distinguished basis of states, which minimalizes quantum effects. That is, the 
quasiparticle creation and annihilation operators in this basis have an appearance devoid of 
polarization clouds. For the XXX and XX Z models with underlying group si (2) the bases 
in question were found by Maillet and Sanchez de Santos |IJ] through the construction of a 
generalized Drinfeld twist ||. 

The ensuing representation of the quantum monodromy matrices coincides, as noted by Terras 
PI, for the case of the rational XXX model with the representation provided by Sklyanin's 
functional Bethe ansatz method ||. An obvious generalization of Sklyanin's method (substi- 
tuting polynomials in the spectral parameter by polynomials in the exponential of the spectral 
parameter) leads us to the conclusion that an analogous coincidence holds true for the trigono- 
metric XX Z model. 

The purpose of the present communication is to report on the generalization of the above 
results to the elliptic sl(2)-XYZ model. For this sake we will make use of Baxter's map of 
the XYZ model onto an ice type model || (which is akin to a cyclic solid on solid model 
H). This brings us formally near to the XXX and XX Z models and allows us to use the 
technique developed in for the handling of the analogous problem in the sl(n)-XXX model. 
The plan of the paper is as follows : Section 2 provides a short survey of the XYZ model and 
its reformulation as an ice-type model. Section 3 deals with the factorizing twists and the 
computation of the operator valued entries of the monodromy matrix. Section 4 contains the 
conclusions. The appendix is devoted to proving some relation needed in Sect. 3. 



2 XYZ model and its relation to ice-type models 

In the framework of the Algebraic Bethe Ansatz M the XYZ-model is determined by the 
elliptic solution of the Yang-Baxter equation 

-^12(^1 — A 2 )-Ri3(Ai — A 3 )i?23(A2 — A 3 ) = R2z{^2 — As)-Ri3(Ai — A 3 )-Ri 2 (Ai — A 2 ) (1) 
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with the notation H{u) = $1 g) , 6(m) = $4 g) and $4 = ^2 
—iq4,e l ' KZ 9/±(z + r/2, g) are the standard theta-functions of a single complex variable.[] 
The somewhat different parametrization as compared to [0] is due to the normalization in order 
to achieve unitarity for the R-matrix. 



The monodromy matrix T(A, {Aj}) (generalized to the inhomogeneous chain ||, [K|) is given 
as the ordered product of Lax operators L^(A — Aj) = Roi(X — Aj) 

T(A, {A,}) = M A -,„)... L 2 (A - z 2 )L l( A - „) = ( ^ {$> {$> ) • (4) 

The presence of the Boltzmann weight d in Eq. @) reflects the non-conservation of spin, which 
is responsible for the absence of a local vacuum for the Lax operator associated with the above 
R-matrix. 

To circumvent the problems arising from the eight vertex nature, we use the vertex-face map 
established by Baxter || to obtain a XX Z type (six vertex) R-matrix by exploiting the relation 

R XyZ (X - (JL)<j) l)V (g> Z m t jV = ^ W ( m ' m '\ l > 0<Am,m' ® Z mt i (5) 

m 

valid for all integers /, V, m,m' such that |Z — l'\ — \m — m'\ = 1 and the summation on the 

r.h.s. is over integers s.t. \m — m'\ = 1, \m — l\ — 1. 

1 For a concise introduction into the theory of theta and elliptic functions we refer the reader to the appendix 
of [pi, whose conventions we shall use throughout our work. 
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The two dimensional vectors 0, z are given by 



M+ i = X(si + n); zi +lt i = X(si + A) 

= X(ti +1 - /j,); zi-x,i = X(ti - A) (6) 

with X(u) = (q^a) and the abbrevation si = s + 2rjl, where s,t are arbitrary complex param- 
eters. 

Relation (||) provides a change of basis in which the Lax operator has a local vacuum indepen- 
dent of the spectral parameter A. There exist a whole family of bases labeled by the integer Z 
and the parameters s, t, which achieve this goal. 
The six weights w(m,m'\l,l') are represented graphically in Figure 1. 
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FIG. 1 The six resulting Boltzmann weights. 

The Boltzmann weights are subject to the star-triangle relation 

^ w(l 2 , h\h, m)w(l 3 , Z 4 |m, l s )w(m, l 5 \h, Z 6 ) = ^ w(Z 3 , Z 4 |Z 2 , m)w(l 2 , m|Z x , l 6 )w(m, Z 4 |Z 6 , Z 5 ) . (7) 



They are parametrized by (h{u) = S(0)H(u)Q(u); ui = + 2r}l - K))\ 

a', = 1; 



a, 



h(\)h(ui- 
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These weights can be arranged into a matrix 



h(\ + 2r])h(uj l ) 

fr(27/)h(m + A) 
h(X + 2r))h(uji) 
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(9) 



which fulfills the modified Yang-Baxter equation JTT| , [fL2] 

i2ia(/ - a 3 )R n {l)R 23 (l - ai) = R 2 3{l)Ri 3 (l - ^ifotf) ■ 
The monodromy matrix related to this modified Yang-baxter equation is f 

^O,l..Jv(0 = ^W(^ — CTi — . . . — (Tjv-l) • • • -Ro2^ — o"i)-Roi(0 



(10) 



(11) 



where denotes the horizontal auxiliary space (with the asssociated spectral parameter Ao), 
1,...,N label the vertical quantum spaces which span the physical Hilbertspace Hn (with 
associated local inhomogeneities {Ai}), and crj equals ±1 depending on whether the arrow in 
the i-th space is up or down (right/left for the horizontal space). It also sets our convention to 
associate the integer in the right lower corner of the graphical representation with the operator 
(cf. Figure 2). 
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FIG. 2 Elements of the monodromy matrix 

2 We use the convention -3Cq i..jv to denote an operator X represented in space with entries in the tensor- 
product of spaces 1, . . . , N. 
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From follows the equation for the monodromy matrices 

#00' (/ - <T\ - ■■■ - cr N )T (l) ® T ,(l - (T ) = T /(Z) ® T (Z - a ,)R o>(l) 



(12) 



and can be represented graphically as in Figure 3. 
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FIG. 3 The Yang-Baxter equation for the monodromy matrix 



The unitarity relation R 2 \Ri2 = I can be checked directly by using the following relation 

h(u + x)h{u — x)h(v + y)h(v — y) — h{u + y)h(u — y)h(v + x)h(v — x) 

— h(x — y)h(u — y)h(v + x)h{v + x) (13) 

which is a consequence of the equality H(u, k)Q(u, k) = H((l + k)u, 2\fkj{\ + k)) with k the 
nome of the elliptic function (Eq. 15.10.20 of [13|]), and the analogue of (|I3|) for H(u) (Eq. 
15.3.10 of 



The representation of eigenvalues of the transfer matrix (H) through these of (|TT[) is explained 
in reference [§] (One has however to keep in mind that our model (11) differs from that of |§ 
by an additional change of basis in the quantum space). 

We will concentrate in what follows on the computation of a factorizing F-matrix for the 
monodromy matrix (|ll|) for an arbitrary fixed value of I. 
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3 The F basis 



The factorizing F-matrix for two sites denned by the relation F 2 iRn — F 12 is 
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(14) 



[12] 



The proof of the factorization property amounts to checking the same relations as those in the 
proof of the unitarity of the R-matrix above. 

The factorizing F-matrix for N sites (N quantum spaces) turns out to be given by formally 
the same expression as found in |7] for the XXX modelQ 



N 

Pa = i[P? 



Zn) 



(15) 



8=1 



where projects on the ctj-th component in the i-th space and the permutation a a is uniquely 
determined through the conditions 



> a <r a {i) ^ a a(i + 1) > <J a {l) 

aa Q (i+i) > a aa{i) if a a (i + 1) < cr a (i) 



(16) 



An algorithm for finding a a for a given a is described in the appendix. 

The modification of the Yang-Baxter equation (0) enforces a particular rule for the handling 
of the integer valued parameter I in the formation of the intertwining matrix R a (l) (related to 
the permutation a), which can be read off from the modified composition law 



R^{1) = R^ {l+l) {k)R°{l) 

k — I — c<j(i) — ... — cr<x(i-i) 

where <jj is the transposition of i, i + 1, and a an arbitrary permutation. 
R a {l) has the intertwining property 

R a (l)To,i...N(l) = T ^(i)... a ( N )(l)R u (l — (x ) • 



(17) 



3 We thank Frank Gohmann for communicating this representation of Fi.„n, which is slightly simpler than 
that quoted in M. 
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The matrix -F\..jv(7) satisfies the factorizing equation 

Rl.. N (l) = F; ( 1.n)(1)Fi...n(1) ■ (18) 

A proof of the latter equation can be found in [0] - the modification of the composition law going 
along with the parameter I being rather immaterial. An alternative argument emphasizing the 
geometric meaning of fll8|) (loosely speaking a pure gauge field representation of R a on the 
symmetric group Sjy as base space) may be sketched as follows: 

It is sufficient to prove Eq. fll8|) for arbitrary elementary transpositions o = af, Oi : (i, i + 1) — > 

(i + l,i), 

R l l ...N = K^1...N) F ^ ■ ( 19 ) 

A special solution F l of flI9| ) is the two site matrix Eq. (0), embedded into the N-fold 

tensor product 

F i = 

and the most general solution of ( fL9|) is given by 

F u . N = X(l,...,{i,i + l},...,N)F i (20) 

where X(l, . . . , {i, i + 1}, . . . , iV) denotes a non-singular matrix in the quantum space ® N C 2 
which is symmetric with respect to an exchange of the space labels i and i + 1 and of the 
respective local inhomogeneity parameters. So we note that the proof of Eq. ([Tj|) is complete 
if we can show that 

Fi...n (P')' 1 = X[l N) (21) 

is symmetric in the labels % and For this purpose it is convenient to consider the 

permutation group element a a of Eq. ([H]) as being decomposed into products of elementary 
transpositions, which also leads to a factorization of R° a into a product of elementary (two 
site) matrices R ax , and to subdivide the sum ^ Q in Eq. (15) into terms with the transposition 



<7j and into terms without. It is furthermore possible to arrange the order of the transpositions 
such that R ai does occur at most once and then on the very right end of the respective product 
of matrices R° x . Schematically we obtain 

a a' a" 

Each factor P a iR aa ' in the first sum has a counterpart P a iiR cr °>" in the second sum related by 
an interchange of the labels i and i + 
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Keeping this in mind one verifies straightforwardly by inspection that the multiplication of Eq. 
(Hby 

from the right side produces an expression which is symmetric in i and i + 1. 



The operators of the monodromy matrix (|TT| ) in the F basis are obtained by using a recursion 
relation, which enables one to express the monodromy matrix for N sites in terms of that for 

N — 1 sites. Starting from the definition Tb,i..-Jv(0 = ^i...jv(0^o,i...jv(0-^i _1 jv(^ ~~ °o) we use the 
identity F^Z) = F^Z - ±±p)F 1)2 ..jv(Z) with F 1>2 ... N (l) = {Pi + P?T 1>23 ... N ) (for the proof 
of this identity we refer the reader to the appendix) 

f 0> i... N (l) = Fi...n(1) t o,i...n(1)K 1 .n(1 ~ CT o) 
= F 2 ... N (l - l±^)F li2 ... JV (Z)T 0)1 ... A r(Z) J P 1 - 1 .. iV (Z - a )F-^(/ - a - . 



Splitting To,i2...jv(Z) into To j2 ..jv(Z — 0"i)-Roi(Z) we arrive at 
Tb,i...jv(Z) 

= F 2 ... N (l - i±^l) j p li2 ... iV (Z) ^(Z - (7 1 ) J P 2 ... iV (Z - oj) T 0i2 ...^(/ - ajRnil) x 

X #T.J\rO - ff - Ol)-F2...Jv(Z - (TO - F^ 2 N (l ~ "o)F 2 T!v( / ~ °"0 ~ 1 "t ~ ' 



i ? i,2...iv(Z)To i2 ... J v(Z - (7 1 ) J R 01 (Z)Fr 2 1 .. JV (Z - (To) 



Inserting the explicit form of 

*wq = ^.jv (/ - ^) Fi,2...v(z)F 2 :^(/ - ^i) = ( ^/ (0 ^ ° iv(/) ) ^ 

we obtain the final expression 

t °' 1 - iv(0 = ( ClL(0 ^..1(0 ) m " ( J " -o) DiJl - a ) , ^ 

(23) 



This relation can be solved recursively starting with the one site monodromy matrix which 
coincides with the Lax operator = R 0i (||). 



We demonstrate the derivation and possible difficulties and their resolution in the case of the 
operator D 1 ^ N (l): 

( Dl_ N (l - l)b' 01 (l) i \ 

Di - n{ ] ~ [ x Di_ N {i)bi__ N {i + 1) (pi_ N ) 1 (Z + 1) ) ^ 

where the lower left element X is given by 

C\..Al)Dl.Al - l)&oi(0 + Dl_ N {l)Cl.. N {l + l)on(i) - Dl M (l)Cl.. N (l) 

This combination vanishes identically by virtue of the Yang-Baxter equation or in graph- 
ical representation: 
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We thus end with a diagonal matrix of the form 
which can now be solved iteratively 
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The computation of the operators -Bi...jv(0> Ci...n{1) proceeds along the same lines to give 

Bl. N (l - l)b> m (l) 



5i..jv(0 



[i] 



with 



Y = Cl_ N {l)B^ N {l-l)b' S^ (dI.X 1 (l+l)Cl M (l+l) 



and for C 

Ci... N {l) 



z di. n (i - 1) (Di. N y\i - l)c' 01 (l) 
o ci_ N (i) 



[1] 

with 

z = cl. N (i - 1) - &Ln(i - 1) (^ 2 °..v) _1 (i - i)cl.. N (i - 1)4(0 • 

The solution of the recursion relation finally yields the result for the operators of the mon- 
odromy matrix ^we use a slight change in notation: 6(A) = h ^l^ and denote k = I — 



~ h(u l+1 ) N fb(Xo-Xi) 

n,(q; 1 , fc+i-jy ) \ u i 



A(Ao) - -^-^^(Ao-A,)^ ®^ Q 6 -i (A ._ A .) J 



D1 



JV 



C,(Ao) = > cRAo-W 



6(A - Xj)b~ l (Xi - Xj) 

(24) 



i=i \ 7 bl 



The above mentioned basis transformation (|^) amounts to splitting the model into sectors 
with a fixed number of turned spins. To obtain the spectrum of the XYZ model one uses the 
operators Bi(X), Ci(X) to construct eigenvectors in the form proposed by |§ (there denoted by 
B k>l (X) etc.). 
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The remaining operator Ai(Xq) is obtained from the quantum determinant det q T(\o) ||14j| , fl5 



det g T(X ) = M^D fc XAo)A.+i,m(Ao + 2r / )-%^ J B M (A )C' fc _ 1 , m (Ao + 2r / ) 

i ^h(X -X i + 2r 1 ) {Zb) 

As the operator D\ is diagonal it can be inverted easily, and the operator Ci possesses a global 
vacuum which is annihilated by its action. This enables one to compute A\ 

h(uo k+i-N ) I / 1 n 
A(X ) = ^ > - N 1 1 



h(u k ) \^ =1 \0 6(A i -A )-V W 



— Cfc-^Ao-AOcKAo-AO/O 0\ o /SfeM 



6(A - Ai) VO lA/^ V ■ 6(A* - A;) -1 



6(A;-A f ) 



4 Conclusion 



The form of the F-matrix, Eq. ( |15|) and the appearance of the monodromy matrix in the basis 
supplied by the F-matrix, Eq's (|4|), (p6|) , are completely analogous to what has been found in 
[0 and [0 for the rational and trigonometric models. The concrete expressions for A, B, C, D 
are in particular manifestly symmetric with respect to exchanges of the local inhomogeneity 
parameters Aj. The quasiparticle operators B and C are free from polarization effects due to 
non-local exchange terms. 

The argument used in - borrowed from [|lj - concerning the identification of operators cor- 
responding to different entries of the monodromy matrix relied on the sl(n) symmetry of the 
rational model. It is not available for the trigonometric and elliptic model. The recursive 
procedure followed instead in the preceding section is equally applicable to the XXX, XX Z 
and XYZ model. 

It seems rather plausible in view of the formal similarities of the rational, trigonometric and 
elliptic models that some version of Sklyanin's functional Bethe ansatz should also be feasible 
in the latter case as has already been achieved for the XYZ Gaudin magnet in [T0| . 
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Appendix 



In this appendix we derive the identity used in the recursion relation (23) using the factorizing 
F-matrix for N sites (for the sake of simplicity we omit the explicit dependence on I, which 
can be restored by using Eq. flTTl) ): 

Fi...n — N {Zi, . . . , Zn) 

N 

Pa = i[P? (27) 



i=l 



and Pf* projects on the ojj-th component in the i-th space. 

With a given a the permutation a a £ Sn is uniquely determined by the following algorithm: 
Let k be the number of l's in the sequence a%, . . . , aN- For i < k we put cr a (i) = I, where I is 
the label of the i-th 1 in the above mentioned sequence. 

For % > k we put a a (i) = k + I', where I' is the label of the i — k-th 2 in the sequence. 
The permutation constructed this way is the only one satisfying the constraints 

aW(i+i) > «<r Q (i) if a a (i + 1) > a a (i) 

a<7„(i+i) > a<r a (i) if (? a (i+ I) < a a (i) . (28) 

Now we will prove that 

F\2...N = F23...N 

If a £ %2 h as fh e form (1, a 2 , ■ ■ ■ , Qfjy) = (1, a) ; a = (a 2 , . . . , ajv) £ %2~ l > we have P a = P/Ps 
and cr Q = cr^ (here and in what follows we identify the symmetric group Sn-i acting on the 
elements (2, . . . , N) as the subgroup of Sn acting on the elements (1,2,..., N)). Thus the part 
of the sum of ( ^7f) corresponding to such a's exactly gives PlF 2 ...N- 

Now consider the case when a = (2, a 2 , ■ ■ ■ , ccjv) = (2,S). It is easy to check that P a = 
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P\P a and a a = a a a l § . . . o x ^-j- where is the number of l's in a and <Jij are elementary 
transpositions and i = u a {i). We have 

PiP a R aa T lt 23...N = P a Ti, aSt (2)... CTa {N)R aa 

= P a Rl t(7& (k+l) ■ ■ ■ R\,cr & {2)R aa 

= P a R a&a ^- a ^ = P a R aa (29) 
where we have taken into account that a^rpj =... = &$ = 2 and that P^PfRu = PfPfl. 
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